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The large luminosity and the increased kinematical domain accessible at the LHC will allow 
to isolate large numbers of events with multiple parton collisions. Interestingly, the hadron 
is probed in different points simultaneously in the process, which allows to obtain informa- 
tion on the correlations between partons in the hadron structure. The whole experimentally 
accessible information on multiple parton collisions consists in the probability distribution of 
the elementary interactions, while the inclusive cross sections, usually considered in largest 
processes, acquire a precise statistical meaning as moments of the multiplicity distribution. 
Conversely different and more "exclusive" cross sections become experimentally viable, pro- 
viding complementary information on the interaction dynamics. The matter is discussed in 
details in the case of hadron-hadron collisions and the role of parton correlations is outlined 
both in the inclusive and in the "exclusive" cross sections. 
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INTRODUCTION 



in 
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The overall features of the minimum bias and of the underlying event in hard hadronic 
interactions cannot be reproduced without implementing multiple parton interactions (MPI) 
in the event generator [l| On the other hand, the direct measurement of MPI requires 

identifying the final fragments of the multiple process, to reconstruct the kinematics at the 
parton level. The sizable reduction of statistics, caused by the request of large momentum 
exchange in each hard interaction, has restricted considerably up to now the possibilities of 
a direct study of the phenomenon and measurements have been possible only in very limited 
phase space intervalsQ] 0| 0] @] - The large energy in the pp cm. system may however justify 
the expectation that, at the Large Hadron Collider (LHC), MPI will not be important only to 
describe the properties of the minimum bias and of the underlying event, giving a significant 
contribution also, as an example, to the cross section to produce jets with large pt- 

A preliminary requirement for the study of MPI in large pt jet production is the separation of 
the contribution of multiple radiation, as the two mechanisms can produce the same final state. 
When looking at the properties of the minimum bias and of the underlying event the separation 
between the two sources is rather ambiguous. In the case of direct observation of the large pt 
scattered partons, one may however expect that the separation of the contribution of MPI form 
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the leading single scattering perturbative Quantum Cromo Dynamics (QCD) mechanism might 
be done more confidently. On one side hard radiation induces an azimuthal correlation between 
the produced jets, which is absent in the case of jets produced through MPI, on the other the two 
mechanisms have a very different dependence on the lower value of the transverse momentum 
introduced to observe the jets and on the cm. energy, being the MPI a power correction to the 
leading QCD single scattering mechanism. Azimuthal correlations were in fact successfully used 
by CDF to isolate and measure double parton collisions at Fermilab and the statement made 
by the experiment, that some 17% of their sample of double scattering events was contaminated 
by triple scatterings^, justifies the hope that the larger energy available at the LHC will allow 
to isolate events with triple and perhaps also with quadruple parton interactions without much 
ambiguity. 

One should nevertheless point out that, while in the case of double collisions the separation 
of hard radiation and of the MPI contribution may be relatively simple, the best strategy to 
isolate the contribution of multiple collisions, when the number of scatterings is large, may not 
be readily identified. Although it might still be relatively simple, at least for triple scatterings, 
in the case of particular channels, like W production plus jets or Drell-Yan plus jets, with the 
jets in particular kinematical configurations. In more general cases, the strategy to separate 
the contribution of MPI from the contribution of the leading single QCD production mechanism 
deserves a dedicated study and will not be addressed in the present paper. Rather we will assume 
that the increased parton luminosity and the sizably higher value of the momentum exchanged in 
a typical partonic interaction at the LHC will allow the direct identification of a sizable number 
of events, not only with double but also with triple and perhaps quadruple parton interactions. 

Interestingly, one may expect that a detailed study of the phenomenon will allow to ac- 
cess several features of the hadron structure still unexplored and most likely unexpected. The 
anomalously small value of the scale factor, characterizing the measurement of double parton 
interactions by CDF at Tevatron[7]|;8i][9], and the need to use a parton density much different 
as compared with the hadron form factor, in the actual Montecarlo generators[l|0|, are indica- 
tions that unpredicted features of the hadronic structure may play a non secondary role in the 
phenomenon. 

A general frame to identify and organize the novel information in a systematic way is hence 
highly desirable. The present paper represents an attempt to provide such a frame. One should 
stress that our goal is by no means to provide a scheme for a general description of the overall 
inelastic event, so our aim is very different in comparison with Montecarlo models. Our final 
states consists only of the large pt partons generated by MPI, where we assume that one will be 
able to recognize and count the number of interactions in some phase space windows. In addition, 
we are interested only in a small fraction of all inelastic events and only in few particular final 
state observables (like for example the jet's momenta). Our aim is to show how such quantities 
are linked to the correlations of the partons in the hadron structure, while remaining as general 
as possible. So we will not attempt to compare various ideas about correlations, which have 
been explicitly implemented, in particular in Montecarlo models. Rather we will try to identify 
observables which have a direct link with correlations and to establish sum rules which, allowing 
to control the consistency of the scheme at a given order in the number of MPI, will help to 
disentangle the effects of correlations from the effects of unitarization. 

A preliminary comment concerns the hard cross sections accessible in the regime of MPI. As 
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it will be discussed in the next sections, in the case of hadron-hadron interactions the integrated 
inclusive cross sections of the QCD-parton model are essentially the moments of the distribution 
in the number of collisions. The whole information of a probabilistic distribution is either given 
by the set of the different terms of the distribution or, alternatively, by the set of its moments. 
One may hence consider two different sets of cross sections, the inclusive cross sections, given by 
the moments of the distribution, and the "exclusive" cross sections, given by the probabilities of 
the different terms of the distribution in the number of collisions. 

The two different sets of cross sections correspond to two different criteria of selection of 
events. In the case of the inclusive cross section of a If -parton scattering, one needs to collect all 
events where there are at least -fT-parton collisions in a given final state phase space interval. All 
events with iV-parton collisions, where N > K, hence contribute with a multiplicity factor (^) . 
On the other hand the "exclusive" cross section of a if-parton scattering process is measured 
by selecting all events with only if-parton collisions in the phase space interval of interest. The 
cross sections called here "exclusive" are in fact partially inclusive cross sections, since one sums 
over all partons outside the given phase space interval and on the soft fragments. 

By measuring the inclusive cross sections, which are related to the moments of the distribution, 
one has the advantage of approaching the problem in a way ordered in complexity. In the case 
of the "exclusive" cross sections, each term is on the contrary related to the whole set of the 
multiple scattering series, since each term contains the probability of no interaction, which carries 
information on the whole series of multiple collisions. In other words each "exclusive" cross 
section requires, at least in principle, the knowledge of an infinite set of non perturbative inputs 
to be evaluated. On the other hand the "exclusive" cross sections are well defined experimental 
quantities, characterizing the regime of hard interactions. It is hence worth trying to discuss also 
the "exclusive" cross sections in the framework of perturbation theory. 

The simplest case of "exclusive" cross section is the hard cross section itself, defined as the 
contributions to the inelastic cross section of all events with at least one hard partonic interaction. 
In the case of hadron-hadron collisions all terms with correlations can be re-summed, when only 
two-body parton correlations are taken into account 0]. As a result of resummation the effects 
of the correlation terms is rather indirect in the hard cross section, in such a way that the 
information on parton correlations cannot be disentangled in a simple way. Rather obviously a 
similar property holds for all "exclusive" cross sections. 

Isolating and organizing the correlation terms in order of growing importance is nevertheless 
possible also in the case of the "exclusive" cross sections. As a consequence of incoherence 
between the different partonic interactions, MPI are in fact described by a (different) probabilistic 
distribution in each phase space interval and one may reduce or increase the average number of 
partonic collisions by restricting or enlarging the interval in transverse momenta and in rapidity 
where final states are observed. It is hence always possible to give an operational meaning to an 
expansion in the number of collisions. The expansion of the "exclusive" cross sections in number 
of elementary collisions allows to isolate the different correlation terms and to obtain interesting 
consistency relations between inclusive and "exclusive" cross sections. 

The paper is organized as follows. After recalling the main features of the simplest Poissonian 
model, which represents the basis of most of the practical implementations of multiple parton 
collisions, we will generalize the approach by including parton correlations in the picture of 
the interaction through a functional approach. Keeping only two-body parton correlations into 
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account, explicit expressions will hence be obtained both for inclusive and "exclusive" cross 
sections, in the latter case by expanding in the number of parton collisions. A simple and 
experimentally testable connection between the two quantities, at the lowest non trivial order in 
the number of collisions, will be finally derived. 



THE SIMPLEST POISSONIAN MODEL 



A standard way to introduce MPI in Montecarlo codes is to assume a Poissonian distribution of 
multiple parton collisions, with average number depending on the value of the impact parameter. 
The approach was introduced long ago[10][l|[ll| and critically re-discussed recently[]3]. Of course 
the Poissonian is only a starting point for the Montecarlo generators, which are in fact much 



more elaborated. Both the JIMMY 
newest variant of PYTHIA 



13(| and PYTHIA [l| include momentum correlations and the 
14l | imposes also flavor conservation. In the present section we discuss 
the simplest option, namely the case where no correlation is taken into account in a hadronic 
interaction at a given value of the impact parameter. Correlations will be introduced in a rather 
general way in the following section. 

One may start by introducing the three dimensional parton density D(x,b), namely the 
average number of partons with a given momentum fraction x and with transverse coordinate b 
(the dependence on flavor and on the resolution of the process is understood) and one makes the 
simplifying assumption D(x,b) = G(x)f(b), with G(x) the usual parton distribution function 
and f(b) normalized to one. The inclusive cross section for large pt parton production as may 
hence be expressed as 

u s = I G(x)a(x,x')G(x')dxdx' = [ G(x)f(b)a(x,x')G(x')f(b - (3)d 2 bd 2 f3dxdx' (1) 

where pf is a cutoff introduced to distinguish hard and soft parton collisions and (5 the hadronic 
impact parameter. The expression allows a simple geometrical interpretation. Given the large 
momentum exchange, which localizes the partonic interaction inside the overlap volume of the 
two hadrons, one may identify with b and b — (3 the transverse coordinates of the two colliding 
partons while (5 is the impact parameter of the hadronic collision. 

Neglecting all correlations in the multi-parton distributions, the inclusive cross section for a 
double parton scattering an is analogously expressed by 
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G(xi)/(6i)<j(xi,x , 1 )G(x' 1 )/(6i - P)d 2 b 1 dx 1 dx[ x 

xG(x 2 )f(b 2 )a(x 2 ,x' 2 )G(x' 2 )f(b 2 - (3)d 2 b 2 dx 2 dx' 2 d 2 (3 
G{x)f{b)d{x,x')G(x')f{b- (3)d 2 bdxdx'^ 2 d 2 f3 (2) 



2' 

which allows an analogous geometrical interpretation. The expression may be readily generalized 
to the case of the inclusive cross section for A-parton scatterings ajy: 

(JN = j j^{f pc G(x)f(b)a(x, x')G(x')f(b - f3)d 2 bdxdx') * d 2 [3 (3) 
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The cross sections are divergent for p\ — > 0. The unitarity problem is solved by normalizing 
the integrand which, being dimensionless, may be understood as the probability to have a iVth 
parton collision process in a inelastic event. One may hence introduce Pn((3), the probability of 
having N parton collisions in a hadronic interaction at impact parameter (3: 



P N {(3) = e-^Cfl j G(x)f(b)a(x,x')G(x>)f(b-f3)d 2 bdxdx'^a s F(P) (4) 

By summing all probabilities one obtains the hard cross section ahard, namely the contribution 
to the inelastic cross section due to all events with at least one parton collision with momentum 
transfer greater than the cutoff p^ : 



= £ J P N ((3)d 2 P = E / <P0 ^ SF ^ N e-°sm = I d 2 (3[l - (5) 



^hard 

N=l " N=l ■ 

Notice that Chard is finite in the infrared limit, which allows to express the inelastic cross section 
as (Jinei = &soft + ^hard with o so ft the soft contribution, the two terms a S oft and ahard being 
defined through the cutoff in the momentum exchanged at parton level, p\. 

As mentioned the inclusive cross sections are related to the moments of the distribution. The 
property is immediate in the case of the Poissonian model. Working out the average number of 
collisions one in fact obtains: 

(N)a hard = f d^ £ NP N ((3) = [d 2 (3 f; "famf e -*sFW = f <p PasF{(3) = as (6 ) 

J N=l J N=l ' J 

which is the single scattering inclusive cross section, while, more in general, one may write: 



(N(N - 1) . (IV -K + l)) = r^ N^-V ( N-K + l ] pN(m 

J N=l 

= jd 2 [3-^[a s F(P)] K = a K (7) 

A few features of the model are worth to be pointed out. 

- The cross sections of the two sets, inclusive and "exclusive" (the latter being the contributions 
to 

Chard due to the different numbers of collisions) are expressed fully explicitly in terms of the 
same quantity, the average number of collisions at a fixed impact parameter crgF(f3). In a more 
general case, when correlations are taken into account, the two sets of cross sections can provide 
a complementary information on the multi parton structure of the hadron. 

- All inclusive cross sections are divergent in the infrared limit. The relation with the multiplicity 
of collisions shows that the divergence is due to the moments. In other words the cause of the 
divergence is the number of collisions, which become very large at low pt. 

- While all inclusive cross sections become increasingly large at low pt, due to the increasingly 
large number of parton interactions, all "exclusive" cross section where the number of partic- 
ipating hard partons is kept fixed become, on the contrary, smaller and smaller for the same 
reason. 
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- In the model one assumes that the transverse hadron size is the same in each interaction, so 
one is identifying the hadron with its average configuration. Hadronic diffraction shows however 
that hadron fluctuations are not a negligible effect also at LHC energies 15] 16]. It was moreover 
pointed out that, most likely, the importance of hadronic fluctuations grows when the number 
of partonic collisions is large 17] 18]. 



3. A FUNCTIONAL APPROACH 



MPI can be discussed in very general terms by a functional formalism 19] 9]. 

Let W n {u\ . . .u n ) be the exclusive n-parton distributions, namely the probabilities to have 
the hadron in a configuration with n partons with coordinates which represent the variables 
(bi,Xi), being b the transverse partonic coordinate and x the corresponding fractional momentum. 
The scale for the distributions is given by the cut off p™ m that defines the separation between 
soft and hard collisions. The distributions are symmetric in the variables U{. 
One may introduce the generating functional: 



z i J } = '^2— / J{u\) ...J(u n )W n {ux...u n )dui...du n . (8) 

The conservation of the probability implies the normalization condition Z[l] = 1. 
The probabilities of the various configurations, namely the exclusive distributions, are the coef- 
ficients of the expansion of Z[J] for J = 0. The coefficients of the expansion of Z[J] for J = f 
give the many body densities, i.e. the inclusive distributions: 



6Z 
5J{u) 



D 2 (U1,U2) 



J=l 



PZ 
5J(ui)5J(u2) 



(9) 



J=i 



Correlations, which describe how much the distribution deviates from a Poissonian, are obtained 
by the expansion of the logarithm of the generating functional, ^[J] = bx2[J], for J = 1: 



F[J]= /d 1 (u)[J(u)-1]^ + ^;^ ! C n (u x ... u n )[J( Ul )-l] ... 

•* n=2 n ' J 

... [J(u n ) - l]dui . . .du Tl 



(10) 



Obviously one has ^[1] = and, in the Poissonian case, C n = 0,n > 2. 

Given the general expressions for the multiparton distributions one can write down an expression 
for the semi-hard cross section which takes into account of all possible MPI: 



@hard 



dp 



Y, 1 

^ n 

n 

£- 

z — J rn 



n\ §J(ui) §J(u n ) 



Z A [J] 



m\ 5J'(u[-P) "' 5J'(u' m - p) 



Z B [J'} 
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n m 

{ 1 - n n t 1 - <m«> «o] } n dndu/ 



i=ij=i 



J=J'=0 



fill 



Here /3 is the impact parameter between the two interacting hadrons A and B and &ij is the 
probability for the parton % (of A) to have an hard interaction with the parton j (of B). 
The cross section is obtained by summing all contributions due to all different hadronic config- 
urations (the sums over n and m). For each pair of values n and m, one has a contribution to 
&hard when at least one hard interaction takes place (whose probability is represented by the 
term in curly brackets in the equation above). 

The cross section, which is analogous to the expression of the inelastic nucleus-nucleus cross 
section in the Glauber modelf^, takes into account both disconnected interactions (which imply 
n = m) and rescatterings (when n ^ m). 

A much simpler expression may be obtained by neglecting all rescatterings. To that purpose one 
may write the term in curly brackets as: 



S = 1 — exp \ ln(l — dij) = 1 — exp 



^2(°a + \ d ii d ij + ••• 



[12) 



rescatterings correspond to repeated indices. Only the first term of the expansion of the logarithm 
hence contributes and all rescatterings are removed by the substitution: 



l-exp^(^) ^J2 &i 3~2^ J2 

ij ij ij k^i,l^j 

The resulting cross section is expressed in a compact way 



Oij<?kl ■ ■ ■ 



(13) 



VhardiP) = exp (d) ■ exp (d') 1 - exp(-«9 • a ■ d') Z A [J]Z B [J'\ 
1 - exp (-d -a-d')] Z A [J]Z B [J'\ 



J=J'=0 



j=j'=i 



(14) 



where all convolutions are understood. In the simplest non-trivial case all correlations C n with 
n > 2 can be neglected and the cross section is more explicitly expressed by: 



PhardiP) 



exp|— J dudu'dja(u,u f )dj>^ 



■ exp 
• exp 



DA(u)J(u)du H — / CA(u,v)J(u)J(v)dudv 



DB(u)J(u)du + - / CB(u,v)J(u)J(v)dudv 



which can be worked out explicitly. One in fact obtains 19] [9] 



(15) 



J=J'=0 



VhardiP) = 1 - exp -- ^a n - ~y^6 ra /n 



(16) 
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where 



and 



a n = (-l) n+1 J D A (u 1 )a(u 1 ,u' 1 )C B (u' 1 ,u 2 )a(u' 2 ,u 2 )C A (u 2 ,U3) . . . 

n 

■■■a(u n ,u' n )D B {u' n )Y[ du i du 'i ( 17 ) 



b n = (-l) n+1 J CA(«n,«l)a(«i,«i)CBK, 

n 

■ ■ ■ C B {u' n _ 1: u' n )a(u' n ,u n ) Y[ duidu'i (18) 

i=i 

Notice that if one works out the moments of the distribution in the number of collisions one 
obtains the same result of the simplest Poissonian model. From Eq.14 one may express the hard 
cross section as a sum of MPI: 



<Thard(P) = [l - exp(-9 • a ■ d')]z A [J]Z B [J' 

= t {d ' a N f )N e- d ^Z A [J]Z B [J'] 



N=l 



J=J'=1 



J=J'=1 



(19) 



The average number of collisions hence is 



(N)a hard ([3) = J2 



N(d-a-d') , & . d , 



N=l 



m 

oo 

= d, h ■ o ■ d j; 

N=0 

= {dj.-a-djJZAlJjZslJi} 



v Z A [J\Z B [J'\ 
N 



J=J'=1 



Z A [J]Z B [J'} 



J=j'=i 



1 



J=j'=i 



D A (xi;bi)a(xix' 1 )D B (x' 1 ;bi — P)dxidx[d b\ = as(P) 



(20) 



where now a{x\x' l ) is the parton-parton cross section integrated with pt > y% since, given the 
localization of the interactions in transverse space, the parton-parton interaction probability has 
been treated as a 5 as a function of the transverse coordinates: a(u,u') = &(x,x')S(h — b'). 
Analogously one obtains 



(N(N- l)) 
2! 



1 



Phardifi) = y J D A (x 1 x 2 -b 1 b 2 )a(x 1 x' l )a(x 2 x 2 ) 

xD B (xix' 2 ; b\ — P,b 2 — j5)dx\dx\d 2 b\dx 2 dx' 2 d 2 b 2 
= <td(P) 



(21) 



and in general 
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(N(N-l)...(N-K + l)) _ 
& hard {P) 

= j D A (xi ...x K ;h... b K )a(x 1 x' l ) . . . u{x K x' K ) 

xDb(x'i . . . x' K ; bi — P . . .bx — (5)dx\dx' x d 2 b\ . . . dxxdx' K d 2 bK 
= <r K (J3) (22) 

which hence shows that when rescatterings are neglected, for any choice of multiparton distribu- 
tions, the inclusive cross sections are given by the moments of the distribution in the number of 
collisions. 



4. "EXCLUSIVE" CROSS SECTIONS 



As discussed in the previous section, the inclusive cross sections are basically the moments 
of the distribution of the number of collisions. The moments contain the whole information on 
the distribution in the number of collisions ordered in complexity: The average number is the 
most basic information on the distribution and is given by the single scattering inclusive cross 
section of the QCD parton model which, precisely for this reason, represents also the simplest 
theoretical quantity to evaluate. Its expression is in fact particularly simple, the convolution of 
the parton distributions (which represent the average number of partons in the hadron in a given 
kinematical configuration) and the single scattering partonic cross section. 

Analogously the ET-parton scattering inclusive cross section gives the iTth moment of the 
distribution in the number of collisions and is related directly to the iiT-body parton distribution 
of the hadron structure. 

As already pointed out, a way alternative to the set of moments, to provide the information 
of the distribution, is represented by the set of the different terms of the probability distribution 
of multiple collisions. Correspondingly, in addition to the set of the inclusive cross sections, 
one may consider the set of the "exclusive" cross sections. As discussed in the previous section 
the case where only two-body parton correlations are taken into account can be worked out in 
details. In the following we will hence work out explicit expression for the single and for the 
double parton scattering "exclusive" cross sections. 

To obtain the distribution of multiple interactions one may start from the partonic interaction 
probability 



n 

where in the product each index assumes a given value only once, in such a way that possible 
re-interactions are not included. The probability of a single interaction may hence be expressed 
as 



d_ 

dg 



n 0--9&ii. 



= ^2°ki n {l-g&ij) 

g=l kl ijy^kl 



(24) 



9=1 
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while the probability for a double collisions is 



1 



a 



2! I dg 



9=1 



2! 



X! X VklOrs J | (1 - gOij) 



kl rs 



ij=£kl,rs 



(25) 



9=1 



The expressions for the single and double "exclusive" integrated cross sections are 



d_ 

dg 



2! I dJ 



9=1 



9=1 



9=1 



i{[X'( 5 )] 2 -X"( 5 )}e-^) 



(26) 



(27) 



9=1 



where X = ^(X] a « + b n /n) with a n and b n given by Eq.'s 17 and 18. 

The physical description of the process is purely probabilistic. As a consequence, the hard 
cross section may be unitarized in each given phase space volume, which hence defines also all 
integration limits. 

By restricting the rapidity window or the pt interval one may control the importance of the 
terms with different powers of o. In a small window only the terms linear in u are relevant. It is 
hence meaningful to expand the terms above in powers of a. Going up to order a 2 one obtains 
for X' and X" 

X'(u,u') = D A (u)a(u,u')D B (u) 

— J Da(u)&(u,u')Cb(u' \u' 1 )a(u' 1 ,ui)DA{ui)du\du' l + A <-> B 

— J CA(u±,u)a(u,u f )CB(u ,Ui)a(ui,ui)duidu± (28) 

X"(ui,u 1 ,U2,u 2 ) = - DA(ui)a(ui,u[)CB(u[,v^)a(u2,u' 2 )DA{u2) + A <-> B 

-C A (u 1 ,U2)a(u 1 ,u' 1 )C B (u' 1 ,u 2 )a(u 2 ,u 2 ) (29) 
At order a 2 , the single scattering "exclusive" cross section d\ is given by 



&i(u,u') = Da{u)(t(u,u')Db{u') 1 — J DA(ui)a(ui,u'i)DB(u'i)duidui 

D A (u)a(u, u')Cb(u', u' 1 )cj(u / 1 , u\)D A^ijduidu^ + A 

Ca(ui,u)&(u, u')Cb(u' , Ui)a(u'i,ui)duidui 
while the expression of the double scattering "exclusive" cross section a 2 is 



B 



(30) 



a 2 (ui,u' 1 ,U2,u' 2 ) = — | D A (ui)a{ui,u\)D B {ui) ■ D A {u 2 )a{u 2 ,u' 2 )D B {u' 2 ) 

+ \p A (u)d(ui,u' 1 )CB{u'i,u' 2 )d(u 2 ,u' 2 )D A (u 2 ) +A^B 
+CA(ui,u 2 )a(ui,u[)CB(u[,u 2 )a(u 2 ,u' 2 )} 



(31) 
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One may immediately verify that 



J ai(u,u)dudu' + 2! J a 2 (ui,u'i,u 2 ,u 2 )duiduidu 2 du 2 = J DA(u)a(u,u')DB(u)dudu 

= (N)a hard = a s (32) 

By restricting the phase space interval of the observed final state, u\ is well expressed by the 
term linear in a, which means that in that limit the probability of interaction is well approximated 
by the average of the distribution, while a 2 is negligibly small. In other words the single parton 
scattering "exclusive" cross section is well represented by the single scattering expression of 
the simple QCD parton model. When the phase space volume is increased, the single parton 
scattering "exclusive" cross section becomes increasingly different from the prediction of the QCD 
parton model and, as it will be shown in the explicit examples of the next two paragraphs, the 
difference allows a direct measure of the importance of correlations. Notice that, by summing 
the expressions of d\ and 02 at order a 2 with the proper multiplicity factors one obtains, as 
shown above, the inclusive single parton scattering cross section, which is correctly given by the 
QCD parton model expression. By comparing the measured inclusive cross section with a sum 
of the measured "exclusive" cross sections, taken with the proper multiplicity factors, up to a 
given order in the number of collisions, one hence has a direct indication of the importance of 
higher order unitarity corrections in a given phase space interval. 

To show how these ideas may be explicitly implemented we work out two explicit examples 
in the following paragraphs 4.1 and 4.2. The transverse parton coordinates are not accessible 
experimentally. On the other hand while, at small x, partons may not be correlated in momentum 
fraction, as conservation constraints may be diluted when the parton population is large, they 
are surely correlated in their transverse coordinates, since they are all confined inside the same 
hadron. A simplest possibility is hence to assume that the only relevant correlations are those in 
the transverse coordinates. To have some indication we will consider two different analytically 
solvable models of the parton densities in transverse space, the Gaussian and the exponential 
model, while neglecting in our two examples all effects of longitudinal correlations. 



4.1. Gaussian model 



One may assume Gaussian distributions for the parton densities and for the correlations: 



D(x,b) 
f(b) 

C(x 1 ,x 2 ;b 1 ,b 2 ) 
h(h,b 2 ) 
h{b,X) 

9(b,R 2 ) 



G(x)f(b) 
g(b,R 2 ) 

G(x 1 )G(x 2 )h(b 1 ,b 2 ) 
c • g{B, R 2 /2)h(b, A) 



1 



exp(-b 2 /R 2 ) 



(33) 



where 
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B = [bi + b 2 ]/2 

b = [bi - b 2 ] (34) 



in such a way that the following relations hold 



J d 2 bg(b,R 2 ) = 1, J d 2 b 2 g(b 1 -b 2 ,R 2 )g(b 2 ,R 2 )=g(b 1 ,R 2 + R 2 ) 

J d 2 bh(b 1 ,b 2 )d 2 b = (35) 

The connection with the root mean square hadron radius (r 2 ) = R 2 is R 2 = 2R 2 /3. We define 
the correlation r c as the value of b where the function h(b, A) changes sign. In this case one 
hence has r c = A. In order to define unambiguously the correlation strength we choose the 
normalization 

f h(b,X 2 )d 2 b = 1 , 

J\b\<r c \ 

since the integral over the whole 6-space gives zero. With this choice rj = e, where e is the Euler 
number. The integrations on the transverse variables of the various terms in Eq.s 30 and 31 give 

D A &D B -> J d 2 bd 2 (3g A (b - frRDgB^R 2 ) = 1 
D A aD B ■ D A aD B -> J d 2 hd 2 b 2 d 2 ! ^(bi - 0, i? 2 )5B(bi, ^)^(b 2 - /?, # 2 )<?B(b 2 , 



2vr(i?f + i^) 

D A aCB&D A -> J d 2 b 1 d 2 b 2 d 2 pg A (b 1 -f3,R 2 )h B (b 1 ,h 2 )gA(h 2 -(3,R 2 l ) 

ce X 2 
~ ix (2R 2 + A 2 ) 2 

CU-rCV _> y d 2 b 1 d 2 b 2 d 2 (3h A (b 1 - /3,b 2 - /3)Mbi,b 2 ) 

_ c^e 2 
~ 4tt A 2 



(36) 



the "exclusive" cross sections <7i and of a 2 are hence explicitly expressed in this case in terms of 
the QCD-parton model single scattering cross section as, the transverse hadron size R and the 
correlations parameters c and A: 
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da i 
dydp t 



da 2 



dyidy 2 dptidp t 2 



da<- 



dydpt 

das x 
dydp t 

1 da s 



<rs 
4vri? 2 

0-5 



1 - 

1 - 

a eff 

da s 



ce 



2as\ z 



tt(2R 2 + A 2 ) 2 



2 2 a S 

c e 



X 



4vri? 2 



+ ce 



2A 2 



2 dyidpa dy 2 dp t 2 

1 do"^ c2<7s 
2a eff dy 1 dp tl dy 2 dp t2 ~ dy x dy 2 dp tl dp t2 



tt(2R 2 + A 2 ) 2 
da D 



4vrA 2 



i 2 2 

+ c e 



4vrA 2 



(37) 



where a e ff is the scale factor which characterizes the inclusive rate of double parton collisions. 
As a function of the root mean hadron radius R, of the correlation length r c and of the correlation 
strength c, the expression is 



1 3 r 3s 2 , , 2 I 

^=8^{ 1 + "' 16 (4T3?F +Ce -3?/ <38> 

where s = r c /R. Notice that neglecting correlation in fractional momenta, in each phase space 
interval a\ is different from a s only by the overall normalization (related to a e ff as shown above). 
At this order in a, a 2 coincides with the expression of the inclusive double parton scattering cross 
section while, in this case, the scale factor a e ff does not depend on x. 

4.2. Exponential model 

A different option is to assume an exponential for the three dimensional parton density 



P(r) = ^-" r , J P(r)d\ = l 
D(x,b) = G(x)f(b) 

f(b) = J p{r)dz = ^bK^b) 

(39) 

where b is the transverse component of r and K\ is the modified Bessel function of the second 
kind. It may be seen that the two distributions Gaussian and exponential are formally connected. 
If one performs a linear combination of the Gaussian distributions (33) with appropriate weight 
one ends up with the expression (39): 



f W eV/4A * I/A)dA ^ 



^/ 4X g(b, l/\)d\ = - V&tfiOife) (40) 



For the next calculations it is convenient to take the Fourier transform of the transverse 
parton density 
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p(b) = wfS e ~ ipt ' h ^ d2pt 



pipt) 



( M 2 +p 2 ) 2 

(41) 



One may assume for the correlation the expression 



C(x 1 ,X2;h,b 2 ) = G(x 1 )G(x 2 )h(b 1 ,b 2 ) 

h(h,b 2 ) = c-f(B)h(b) 

Kb) = j^y J e "iP* • h h( Pt )d 2 Pt 

h(p t ) = v 



d K 4 7 2 



= 2^712^ (42) 



in such a way that 

fr(0) = y /i(6)(i 2 6 = (43) 
As a function of the coordinates the expression is 

h(b) = A-^bK^Kb) - ^K 4 b 2 K 2 ( K b)} (44) 

where K\ and K 2 are modified Bessel functions. The root mean square value correlation term 
is performed in the same way as in of the density is (r 2 ) = R 2 = 12/ fj 2 and, like in the previous 
case, defining the correlation length as the value of b where the function h(b) changes sign, one 
has r c = Xq/k, where xq ~ 2.386 solves the equation 



4x0^1 (x ) =x%K 2 (x ) (45) 

The normalization of the correlation function is performed in the same way as in the Gaussian 
model, it gives now rj = 3.456. The integrations on the transverse variables of the various terms 
in Eq.s 30 and 31 give 
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D A oD B - J d 2 bd 2 f3f A (h - (3)f B (b) = 1 
D a &Db • D A aD B - y d 2 b 1 d 2 b 2 d 2 /3f A (b 1 - /3)/ B (bi)/A(b2 - /?)/ B (b 2 ) 



D A aC B aD A 



7ttR 2 

dH^hd 2 PUibi - P)h B (h 1 ,h 2 )fA(b 2 - (3) 



3r ? Z7 
7TR 2 



x 2 R 2 



12r? 



F(a) ee 



3 + 44a - 36a 2 - 12a 3 + a 4 + 12a(2 + 3a)ln(a) 



3(a- l) 6 a 



d 2 b 1 d 2 b 2 d 2 (3h A (b 1 -/3,b 2 - 0)Mbi,b2) 



2 2 

c ry 



30vr r-2 



Notice that the function F(a) is regular near a = 1, the expansion being 



F(o) 
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(46) 



(47) 



As a function of the root mean square hadron radius R, of the ratio s ee r c /R and of the 
correlation strength c, the expression of the effective cross section is 



{ 1 + c ,. 14 . F (_L£f)+ cV ^(i|)} 



T eff 



7llR 2 



(48) 



One may hence obtain information on the correlations by comparing the measured "exclusive" 
single scattering cross section with the measured inclusive single scattering cross section: 



dai _ da s = da s as ^ 
dydp t dydpt dydp t a e ff 

The cases considered here are relevant for a phase space window where multiple collisions, of 
order higher than two, do not give important contributions; condition which may be always 
achieved by taking, for example, the lower cutoff in p t large enough. Notice that the importance 
of higher order scattering terms, namely of unitarity corrections, is controlled irrespectively of 
any simplifying hypothesis on correlations, by verifying the validity of the sum rule Eq.32. 
The relation above allows to obtain the value of cr e ff- By subtracting from the effective cross 
section the term 3/(77r^ 2 ), or 3/(87i\R 2 ), depending on the prejudice on the parton density in 
transverse space, one obtains a direct information on correlations. Notice that, although the 
two options considered, the Gaussian and the exponential, correspond to rather different parton 
densities, the quantity to be subtracted is basically the same in the two cases and it is essentially 
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determined only by the value of the mean square hadron radius R 2 . Indications on the value of 
R 2 are obtained by the measurement of the generalized parton distributions 22] 23]. In the case 



of gluons one may estimate R 2 ~ .42 fm 2 24] 



One should point out that the validity of Eq.49 is easily tested, since the dependence of the 
difference in the right hand side of the equation, integrated in a given phase space interval, on 
the rapidity window and on the lower cutoff in pt has to be the same as the dependence of <r|, 
integrated in the same phase space interval. So, if the sum rule of Eq.32 is satisfied while the 
relation in Eq.49 is not, one may safely conclude that one has direct evidence of longitudinal 
correlations in the hadron structure. 

The analysis may be repeated by including higher order multiple scattering terms in Eq.30 
and in the following. In such a case, one may write relations similar to Eq.49, linking higher order 
multiple scattering inclusive cross sections with various combinations of "exclusive" cross sec- 
tions. Independent constraints on the correlation parameters may hence be established, hopefully 
allowing an unambiguous identification of the correlation parameters. 



5. CONCLUDING REMARKS 

MPI are originated by the large flux of partons in high energy hadronic collisions. The 
dominant contribution to the process is hence due to the terms which maximize the incoming 
parton flux. Given the number of partonic collisions, the incoming flux is maximized by the 
set of terms where elementary hard interactions are all disconnected one from another or, more 
precisely, where the connection is only due to soft exchanges. Being the interactions uncorrelated 
with respect to hard dynamics, all correlations observable in the final state are originated by the 
initial state flux, namely by the hadron structure. In such a scheme all partonic collisions add 
incoherently in the final cross section, leading to a purely probabilistic picture of the process. It 
is hence natural to identify two different classes of hard cross sections, the inclusive cross section, 
related to the moments of the distribution, and the "exclusive" cross sections, related to the 
different terms of the probability distribution in the number of hard collisions. 

A relevant property is that the picture of the interaction is probabilistic in each phase space 
interval, in such a way that one may associate a different probability distributions to each different 
phase space choice to observe the final state. Including the possibility to select events which are 
fully or only partially contained in the given phase space interval. Any choice identifies an 
interval in pt and restricts in a precise way the limits in x of the initial state, which identifies the 
domain of definition of the probability distribution of multiple collisions. One should point out 
that the same restricted intervals define also the integration range of the integrated terms, which 
appear in the "exclusive" differential cross sections (e.g. the integrations in Eq.30). Integrated 
terms appear in the "exclusive" differential cross sections because of the normalization of the 
probability distribution, their origin may be traced back to the virtual corrections of the process. 
Normalization hence fixes unambiguously the integration limits of the virtual terms, which must 
coincide with the kinematical limits imposed to the real terms by the choice adopted to select 
the final state. 

Another relevant point is that, as proved in the appendix, the components of the hadron 
structure, namely the terms D and C n , are well defined quantities in the picture of the interaction, 
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as they do not mix when the selection choice of the final state is modified. The effect, to modify 
the selection choice of the final state, is in fact only to change the integration domain of each 
term, namely to probe the multi-parton structure of the hadron in different domains in x and 



As it was pointed out, the inclusive cross sections can provide a direct information on the 
correlations of the multi-parton structure of the hadron, while to evaluate the "exclusive" cross 
sections one needs, in principle, an infinite non perturbative information. The expansion of 
the "exclusive" cross sections in the number of parton collisions allows however to identify the 
dominant terms and to isolate the leading contributions due to the presence of correlations. 
Explicit relations between inclusive and exclusive cross sections may hence be written at each 
order of the expansion in the number of collisions, while the importance of higher order terms 
may be controlled by selecting the phase space interval where to observe the final state. The 
comparison between experimentally measured inclusive and "exclusive" cross sections, at a given 
order in the number of collisions, provides a model independent tool to judge on the importance 
of higher order unitarity corrections (namely by checking the validity of sum rules like Eq.32). 

In the present paper we have worked out in details, at the lowest non trivial order, the 
simplest case where only the correlations in the transverse coordinates are taken into account. 
Of course correlations will depend also on fractional momenta and on the different kinds of 
partons considered (gluons and quarks, distinguishing different flavors and the valence from the 
sea). To some extent the different kinds of partons involved may be selected by looking at definite 
final states produced and at the corresponding values of momentum fraction involved (a prompt 
photon plus five jets at low x are likely to be generated by a triple scattering involving a quark 
and five gluons etc.). The information on correlations will hence to be related to the different 
kinds of initial state partons involved in the interaction. 

While the importance of higher order unitarity corrections is tested, in a given phase space 
interval, by verifying the sum rule in Eq.32 (and similar when including a higher number of 
collisions) terms with correlations are isolated by using relations like Eq.49, which allow to obtain 
the value of the effective cross section and of similar scale factors for higher orders unitarity 
corrections. The importance of correlations is tested by comparing with the expectation in the 
uncorrelated case and, in particular, longitudinal momentum correlations will show up with 
the dependence of the effective cross section (and of the similar scale factors for higher orders 
unitarity corrections) on the rapidity window. With respect to longitudinal correlations one 
should anyhow point out that, while indications on the effects of longitudinal correlations may 
be obtained along these lines in pp collisions, a model independent separation of longitudinal 
and transverse correlations may be obtained by looking at MPI in the case of hadron-nucleus 
collisions! 2 11. 



As already mentioned in the main text the "exclusive" cross sections are defined with respect 
to intervals of variation of the kinematical variables, which are to some extent arbitrary. Even 
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though the result may look intuitive, we like to analyze formally what happens when varying the 
intervals of the kinematical variables: Starting with the formalism stated in eq. (8) we call U the 
field of variation of u and divide it into two mutually exclusive parts U = TZLiS. Correspondingly 
it is convenient to denote by r the variable u when if varies within 7Z and s when it varies within 
S. Then for the exclusive distribution within 1Z we have: 

Wi(n) = W 1 (r 1 )+ [ W 2 (n,s)ds+l [ W 3 ( ri ,s,s')ds ds'+^- [ W 4 (ri,s,s' ',s")ds ds' ds" + ... 
JS * Js 3- JS 



W 2 (n,r 2 ) = W 2 (n,r 2 ) + J W 3 (n,r 2 ,s)ds + ^ J W 4 (r 1 ,r 2 ,s,s')ds ds' + . . . 

These relations may be summarized by means of the following procedure: Let us define the 
projectors R(u) = 1, S(u) = for u € 1Z R(u) = 0, S(u) = 1 for u £ <S, R + S = 1; and the 
generating functional Z[J] = Z[R ■ J + S]. Then expanding Z[J] as in eq (8) and using the 
symmetry of the W n it is easy to attain the result: 

^[ J ] = J2 f J ( ri ) • • • J^nWnin . . . r n )d n ...dr n 



Now we perform the same substitution in the logarithm of the generating functional F[J] = 
T[R ■ J + S] and using 1 — S = R the projector R is found in front of every integrand, so the 
overall effect is a restriction in the integration domain and the equivalent of eq. (10) is easily 
found: 



T[J]= [ D 1 (r)[J(r)-l]du + f I / C n (n . . . r„) [J(n) - l] . . . [J(r„) - l] dr x 
Jn £^ 2 n]J K 



, dr ri 



The conclusion is hence that the inclusive quantities are not modified by the variation of the 
interval where the kinematical variables are defined, on the contrary the variation of the exclusive 
quantities can be worked out. With respect to the previous formula the difference arises only 
from the different integration domain of the terms D and C n . 
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